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Abstract 

We determine the wavefunctions of electrons bound to a positively charged mesoscopic 
metallic cluster covered by an insulating surface layer. The radius of the metal core and the 
thickness of the insulating surface layer are of the order of a couple of angstrom. We study 
in particular the electromagnetic decay of externally located electrons into unoccupied 
internally located states which exhibits a resonance behaviour. This resonance structure 
has the consequence that the lifetime of the "mesoscopic atoms" may vary by up to 6 
orders of magnitude depending on the values of the parameters (from sec to years). 

1 Introduction 

Clusters consisting of some 10 2 to 10 5 atoms, which can be produced for instance in the course 
of the condensation of vapour, have attracted great interest during the last 20 years. These 
"mesoscopic" clusters were found to exhibit classical as well as quantum mechanical features. 
In the special case of metallic clusters , the conduction electrons, which move in the average 
potential produced by the positively charged ions, exhibit shell effects. As a consequence, the 
binding energy of the cluster becomes a function of the number of free electrons (and thus also 
of the number of ions). In turn, this implies measurable variations of the formation probability 
of a cluster as a function of its mass in a thermal vapour-gas equilibrium. 

Other interesting types of mesoscopic systems are the so-called "fullerenes" which represent 
a crystal of atoms covering a closed, usually spherical surface [2|. A number of review articles 
exist on this exciting new field of physics (see for instance Ref. |3j and Ref. [I] and books 
mentioned at the end of Ref. |T] and Ref. [2] ) . 

In the present paper, we consider a mesoscopic system which consists of a spherical metallic 
core of radius R2 carrying Z positive charges and surrounded by an insulating layer of external 
radius R\ (see Fig. 1). We determine analytically the bound states of a single electron in 
the potential V(r) which is produced by the positive surface charge density of the metallic 
core and the insulating dielectric surface layer. As in the case of an ordinary atom, these 
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single electron wave functions describe also approximately the case that several electrons are 
bound in this potential. For Z bound electrons, the system has total charge 0. Because of the 
obvious similarity of this system with an atom or ion, we refer to it as a "mesoscopic atom" or 
"mesoscopic ion". 

As one can surmise, there are two different types of bound electronic states: electrons 
localized in the metallic core inside of the insulating layer (region II, Fig. 1) and electrons 
localized in the vacuum outside of the insulating layer (region 0, Fig. 1). The bound states 
localized outside are reminiscent of the bound states of an electron in the Coulomb potential 
of the atomic nucleus. We are mainly interested in these externally localized electronic states 
("class 1-states") and their decay into lower- lying bound states predominantly localized within 
the metallic core ("class 2-states"). 

The paper is organized as follows: 

In Sec. 2, we formulate the single-particle Hamiltonian and determine the eigenfunctions 
and eigenenergies of bound electrons. In doing so, we take into account the polarization of 
the insulating surface layer which is produced by the uniform positive charge distribution on 
the surface £2 of the metal, but we neglect the additional polarization effects produced by the 
externally localized electrons. 

In Sec. 3 we determine these additional polarization effects within the classical Maxwellian 
electrodynamics. They lead to an additional potential acting on the electrons which is of 
importance close to the surface S\ of the insulator. So far, we have not yet studied the effect 
of this additional potential on the externally localized electrons. 

The Sec. 4 is dedicated to a presentation of the results. In Subsec. 4.1, we describe the 
procedure which we used in order to obtain the eigenvalues and eigenfunctions and we present 
typical examples of spectra and wavefunctions. In Subsec. 4.2, we discuss the electromagnetic 
decay of a predominantly externally localized bound electron into a lower-lying internally lo- 
calized state. We show that the transition rate and the corresponding life-time may vary by 
several orders of magnitude depending on the parameters of the system. 

In general, the transition is hindered by the small overlap between externally and internally 
localized states ("class 1-states" and "class 2-states") if the insulating surface layer is sufficiently 
thick. In particular, we discuss the special role in this decay of " class 3-states" which exhibit 
non-negligible amplitudes both in the external region and the internal region II. 

Finally, in Sec. 5, we summarize the results and discuss some related open questions. 



2 Theory 

In the single-particle Hamiltonian for an electron bound to a positively charged and insulated 
metallic core (M = mass of the electron) 

H = f + V = -^A + V(r) , (1) 

the potential V(r) contains the interaction of the electron with the positive surplus charge Ze 
(eo = elementary charge) and with the ionic background ("jellium") of the metal. Furthermore, 
it contains the effects of the insulating surface layer. The radius R2 of the metal and the 
thickness (i?i — R2) of the surface layer are assumed to be of the order of 5 to 20 A(l A= 
1CT 10 m) and the number Z of positive charges of the metal is considered to be 1 to 5. For 
these dimensions of the system we may expect that, on the one hand, the dynamics of the 
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Figure 1: Schematic view of the charged metallic cluster (radius R 2 , surface Sz) covered by an 
insulating layer (outer radius R\, surface Si). States of bound electrons predominantly localized 
in the surrounding vacuum (R > Ri, region 0) are called "class 1-states" in distinction from 
states localized in the metallic core ("class 2-states"). 



electron is fully quantum-mechanical and, on the other hand, we may use macroscopic concepts 
in formulating the potential V(r), at least as a first approximation. So we assume that the 
positive surplus charge is uniformly distributed over the surface S 2 of the metal 

p 2 (r) = p S2 ■ 5(r - R 2 ) = ^5(r - R 2 ) . (2) 

The Coulomb potential produced by this simple charge distribution acting upon an electron 
(charge — e ) is given by 

- e / = 6 (R 2 - r) - ^6 (r - R 2 ) , (3) 

J \r — r\ R 2 r 

where 6q(x) is the Heaviside function 

o (x) = 



1 for x > 
for x < 



Inside of the metal, the conduction electrons feel not only the constant potential produced 
by the positive surface charge but also the attractive interaction with the ionic background 
which we represent by a constant negative potential — V$ (V^ > 0). 

On the other hand, the average energy gap above the Fermi level in the insulator may act 
upon the electron like a repulsive barrier of a height Bq > 0. 

Last not least, there are polarization effects produced by the positive surplus charge and 
also by the charge distribution p e (r) of the externally (i.e. predominantly in region 0) localized 
electrons. The polarization effects are investigated in Sec. 3 on the basis of Maxwell's equations 
for the electric field E and for the displacement field D in the presence of homogeneous dielectric 
media. In the case that we only take into account the polarization effects produced by the 
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positive surplus charge (2), the following modifications of the Coulomb potential (3) are to be 
considered: Within region I, i.e. within the insulating layer, the Coulomb potential —Zel/r is 
to be replaced by 

7..1 7,:l 7, 

'0 



Zel 



Zel 



Zel ■ (ei ~ 1) 



r r ■ s\ Rx ■ E\ 

and within region II, i.e. within the metallic core, the constant potential —Ze\jR2 is to be 
supplanted by 



Zel 
R2S1 



R\E\ 



Zel + Ze i ( 1 1 

R2 V-R2 R\ 



ei-l 
£1 



Within region 0, i.e. in the vacuum surrounding the insulating layer, the value —Zel/r of the 
Coulomb potential (3) is unchanged. 

One should note that the potential in region I assumes the value —Ze\jR\ at the outer 
surface r = R\ of the insulating layer. 

Let us denote the sum of the constant parts of the potentials in region I and II by 



and 



B := Bn 



K? - Zet 



E\Ri 



(ei - 1: 



1 

Ri 



CI 



(4) 



(5) 



,-R2 tt\J v £1 

If Z is equal to 1 or 2 and the dielectric constant E\ exceeds 1 only slightly, and if the radii 
are in the range of values treated in Sec. 4, the 2 nd terms on the right-hand sides of Eqs. (jlj) 
and (jSj) are much smaller than the 1 st terms. So far, it is only these cases which we treated 
numerically. 

Altogether, the potential Vir) felt by the electron has the form 



V(r) = - 



(R2-r) + 



B- 



z4 

E\T 



Zel 



6 {R l - r)9 (r - R 2 ) °-9 (r - R x ) . (6) 



The constants Vn, E\ and Bq are empirical parameters. The potential V(r) is shown in Fig. 2 
for a specific choice of the parameters. 

The eigenfunctions if)(r,'&,(p) of the Schrodinger equation 



2M 



Aip + V(r)ip = Eip 



separate 



4>(r,0,<p) = R(r)Y lm (#,(p) 



u{r) 



and the radial functions R(r) and u(r) have to satisfy the differential equations 



2M \dr 2 
h 2 ( 



id 1(1 + iy 



r dr r ' 
d 2 1(1 + 1) 



R + V(r)R = ER 



2M V dr 2 



u + V(r)u = Eu 



(7) 
(8) 

(9) 
(10) 
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Figure 2: Semi-quantitative picture of the potential V(r) for the following choice of parameters: 
R 2 = 20 A; R x = 30 A; Z = 1; W = 3.3 eV; V u = 6.3 eV; 5 = 3 eV; e x = 3 ; Zeg/(e r • i? 2 ) = 
0.24 eV; Zeg/(ei ■ = 0.16 eV; Ze 2 /Rx = 0.48 eV . 



respectively. 

For r < R 2 (region II), the potential V(r) is constant and thus the Eq. represents the 

Of the 2 linearly independent solutions, 
Up to a normalization factor which 



differential equation for "spherical Bessel functions'' 
only the regular one, ji, qualifies as a physical solution, 
will be determined later on, the radial function R(r) in region II is simply given by 



Rli(r) = ji(Knr) 



where Ku is defined by 



2M 



\ s 2 



\E\ 



111 



(12) 



and ji(x) is the spherical Bessel function of the 1 st kind (see Chapt. 9 of Ref. 5J). As bound 
states have an energy between the well depth — (Vn + Zel/R 2 ) and 0, the parameter Kn is 
assured to be real. 

For i?2 < r < R± (region I) let us postpone the treatment of the general case and first turn 
to the simpler situation when the r-dependent term Ze^j[e\r) is negligible in comparison to 
the constant B in the entire region I. This is so, if the inequality 



11 



« B 



(13) 



holds. For the case shown in Fig. 2, the condition is fulfilled. Generally speaking, the barrier 
heights B range from 2 to 5 eV, and the dielectric constant E\ may have values between slightly 
more than 1 up to more than 100 (Examples: S\ (caoutchouc) = 4; e\ (Ti02-crystal) = 115). 
Let us now investigate the case that the potential (JBJ) can be approximated by 



V(r) 



Vn + ^ 



R2 



Zel 



0o{R2 - r) + B0 o (Ri - r)9 (r - R 2 ) °-9 (r - Rx) . (14) 



The radial Schrodinger equation (jSJ) in region I then assumes again the form of Bessel's differ- 
ential equation. But as the energy E is lower than the barrier height B, the solution in region I 
is given by a linear combination of two independent "modified spherical Bessel functions" (see 
Chapt. 10.2 of Ref. 0). 
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We choose the pair {(— i) l ji(iy), — (zy)|, where ji(iy) and hf\iy) are the spherical 
Bessel-function and spherical Hankel-function of the 1 st kind, resp., with imaginary arguments. 
The dimensionless real variable y is related to the radial coordinate r by 

y = Kl r, (15) 



Kl = ^(B+\E\) . (16) 

The two factors (—i) 1 and — fi 2 are chosen so that the pair of solutions are simply related to 
the modified Bessel function of integer order I l+ i(y) and K l+ i(y) for which a number of rapidly 
converging series expansions exist (see Ref. [3], Chapt. 10.2) 

{(-iYUwl—^m = • ( 17 ) 

The radial function R(r) in region I thus has the general form 



*« = W^F**^ ■ <18) 



For r > Ri (region 0), it is more convenient to deal with the differential Eq. (fTT)|) for u( 



r 



which has the form 



Ti 2 ,„ . h 2 l(l + l) Ze\ , . 

- T u'\r) + ^-—lu--^u = Eu. (19) 



2M v ' 2Mr 2 r 
We have to find a normalizable solution of ([19)1 . It turns out that the only normalizable solution 
of (|19jl has the form 

u(r) = C z l +1 U(a , b , z )e~^ , (20) 
where the dimensionless variable zq is related to r by 

z = 2\r (21) 

with 



A^F^F (22) 



'2M\E 
H 

and where u> = U(ao, bo, zq) is a particular solution Kummer's differential equation 



d 2 w . , du> 
^o-tt + (°o - ^oH a w = . (23) 



The parameters ao and &o in (pH| are related to the parameters in Eq. (JTHJ) by 

ao = ?i±ihL* (24 ) 
6 = 2(1 + 1) (25) 



Ze^ [2M , . 

d » = irVm (26) 



and the solution U(ao, bo, Zq) of (f2l3|) is denned as follows: 
We write bo in the form 

b = l + n, (27) 

which implies 

n = 21 + 1 . (28) 

Then the function U(ao, 1 + n, z ), where n — 1, . . . reads (see Ref. Eq. (13.1.6) or Ref. [Hj, 
Eq. (13) of Chapt. 6.7.1) 

U(a ,n + l,z ) = -ppT r • {M(a ,n+ l,z ) In z + 

n\i (do — n) 

+ E ri^itll ^( Q ° + - V(i + ^) - V(i + n + ^)] } 

i/=0 1" i" -U^- J 

r(ao) is (! - *)« ^ 

In f|29|) . the function ip(x) is the logarithmic derivative of the T-function 

, . , dT(x) . d\iaT(x) , , 

:= -£i/T(«) = (30) 

and (c), is Pochhammer's symbol 

(c),:=c(c+l)...(c + z,-l) = ^Z) . (31) 
Furthermore, M(ao, 6o 5 ^o) is the confluent hypergeometric function 

r ( a o)r(fe + v\ 

We note in passing that the complicated form of the solution (j2HJ) is due to the fact that the 
parameter b in U(a , b , z ) is an integer. For non-integer b , the function U(a , b , z ) depends 
in a simple way on two confluent hypergeometric functions (see Ref. jH], Eq. (13.1.3)). This 
form of U(ao, bo, zq) becomes singular if bo tends to an integer value. The limiting process which 
leads to the form (J2"§)) is discussed in Ref. 0. 

For any value of b , the asymptotic behaviour of U for 9R(zo) ~~ > 00 is given by 



U(a ,b ,z ) -> 4^7 



i + o i 



(33) 



(see Ref. jS], Eq. (13.1.8)). From (}3*3*j) and (|2Ujl it is clear that w(r) is a normalizable function 
irrespectively of the value of the energy E. However, it turns out that the continuity conditions 
can only be met for certain discrete values of E. The radial function Ro{r) in region is thus 
defined by 

Mr) = = Co(2Xr) l U(ao, 21 + 2, 2Xr)e~ Xr . (34) 
zAr 
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Before we proceed to the explicit formulation of the continuity conditions, we return to the 
form of the solution in region I in the general case that the r-dependent term Ze^/ (efr) is not 
neglected. Then the differential equation for u(r) in region I has the form 

' u + [ B Si \ u = Eu . (35) 



2M \dr 2 r 2 J \ s\r 
Eq. (jHSj) resembles Eq. (fTHj) . In fact, with the replacements 



Eq. (0 becomes Eq. (|331) . 
With the ansatz 



Z - - , (36) 
E^E-B , (37) 



u (r) = ^ +1 w(zi)e-^ (3* 



we obtain again Kummer's differential equation for w(z 

d 2 w dw 
dz 2 Zl dz\ 



z i~Tj2 + ( b o ~ z i)— -aiw = 0, (39) 



where bo is again given by (|2l)|) or (|2"7|) . but ai is obtained from a by the replacements (|HEJ), 

2(Z + 1) - ^ 

°I = o ' ( 40 ) 



Ze 2 I 2M 



* = ^Vi^Bi (41) 



and the dimensionless variable zi is related to r by 



\2M-\E-B\ , s 

zi = 2 • W L_ = 2mr . (42) 



The parameter Kj defined in (I16J) agrees with ^2M/h 2 -\E-B\ because E < and 5 > 0. 
The general solution w of Eq. (|3*H|) can be written as a linear combination of the function 
U(ai, 21 + 2, zi) and, as linearly independent 2 nd solution, the confluent hypergeometric function 
M(a 1 ,2l + 2,z 1 ): 

w(zi) = A 1 M{a l , 21 + 2, z T ) + BiU(ai, 21 + 2, z T ) . (43) 

We note that we had to suppress the solution M(a , 21 + 2, z ) in region because it would lead 
to a non-normalizable solution. 

The radial solution R(r) in region I thus assumes the form 

*M = , (44) 

Ri(r) = {2K ir ) l e' Kir ■ [AiM(a h 21 + 2, 2/qr) + BiE7(ai, 21 + 2, 2«ir)] . (45) 

We now turn to the continuity conditions which enable us to determine the amplitudes A\, 
B\ (or Ai, Bi), Co, and the energy eigenvalues E. Once these parameters will be obtained, we 
shall normalize the entire radial function to 1. 
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The radial wave function R(r) and its 1 st derivatives must be continuous at the limits r = R2 
and r = R\ of the regions II and I: 



Ru{R2) — RiyR-i) ■ 
dR u {R 2 ) dRi(R 2 ) 



(IR2 (IR2 
R\{Ri) = Ro(Ri) , 
dRxiRt) dRoiRj 



(46) 

(47) 
(48) 
(49) 



dRi dRi 
We proceed as follows: 

Eqs. (@n|) an d (|4T|) are 2 linear inhomogeneous equations for the amplitudes Aj, Bj (or Aj, B{). 
We solve them for these amplitudes and substitute the result into the equation 



din Ri(Ri) dkiRo{Ri) 



dR! 



dRi 



(50) 



i.e. the ratio of the Eqs. (}4*9~)) and (}4*Sj) . Eq. (jBTIj) represents a transcendental equation for the 
energy. It is solved numerically. The amplitude Cq can subsequently be obtained from ()48|) . 
The resulting wavefunction R(r) is not yet normalized to 1. This can be easily achieved by the 
replacement 

R{r) -> R{r) := M ■ R(r) , 

with 



M=\J drr 2 R 2 (r) 



(51) 
(52) 



Let us present this procedure in more detail: 
We introduce the abbreviations: 



x=K\r 



(2)/ 



-K\r 



X = K\T 



x=K\r 



ff\2 Kl r) = [f{ 



(1) 



y=2n\r 



y l e-zM(a h 2l + 2,y) 



y=2Kjr 



fi 2 \2K ir ) ee [f} 2) (y)] y =2 Kir := \y l e- V iU(a h 2l + 2,y)\ 

Furthermore, we use the prime ' to denote the derivative with respect to the variable given as 
argument: 



?(2), 



(53) 

(54) 

(55) 
(56) 



ff ] '(2^2) 
fl(KlM 



dx 

' df?\y) 
dy 

djijx 
dx 



y=2KiR 2 



- X=K, n R' 2 
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a.s.o. 



From the continuity equations (|4lj|). (jj^j) we find the amplitudes Ai, -Bi and A[, B\ to be given 
by very similar relations 



, Ji{^iiR2)fi 2) '(«i^) - ^/r(«i^)iK«ni22) 

A! — 



e(2). 



Si 



The Wronskian in the denominators is given by (see Ref. [Sj, Eq. (10.2.8)) 

W(/, (1) (aO, // 2) (a;)) := /«(*)// 2) '(*) - /^(s)/^ '(*) = 



r(2). 



(57) 
(58) 

(59) 
(60) 

[W(fT\y)jT(y))}y=2 Kl R 2 ' (61) 

The Wronskian in the denominators of ()60)1 and (|61J) can be related to the Wronskian W{M(a\, 21+ 
2,y),U(a h 2l + 2,y)): 



2x 2 



ji{K ll R 2 )fi 2) \2K l R 2 ) - ^-jr>(2 Kl R 2 )f l (K ll R 2 ) 



(2), 



[w(P(y),fr(y))]y 



2Kj ■ 
(2), 



-2k x R 2 



ji(K ll R 2 )ff ) '(2K 1 R 2 ) - ^-fl L) {2K 1 R 2 )j[{K 11 R 2 ) 



2ft! • 



W(fT(y), ~fl 2 \y)) := /z (1) (y) if '(y) - '(v) 

^(i/%), 7/ 2) (2/) = y 2 'e- y ■ W(M(a r , 2/ + 2,y), U(a h 21 + 2, y)) 
Using the differential Eq. flU)), one finds for W(Af, C/) (see Ref. 0, Eq. (13.1.22)) 

T(2/ + 2) • e v 



(62) 



lU(M(ai,2/ + 2,y),f/(ai,2/ + 2,|/)) 



r(oi) • y 



21+2 



(63) 



and thus for W(f} X \y),f?\y)): 



?(2), 



w(^ {1) (y),/r(y)) 



(2), 



r(2/ + 2) 



(64) 



Substituting (p)9"j) into (|57|) and (|3H)l and, analogously, substituting (jMj) into (|50*|l and (pTj) we 
obtain the following forms of the amplitudes in region I: 



A 1 



7T 



ji(KnR 2 )fi 2) 'foifc) - ^^{^R^k^) 



Bi 



2k\R\ r 



7T 



Hi 



Bi 



" r(2/ + 2) 

4^i?|r(ai) 

r(2/ + 2) 



j l (K ll R2)ff )) '(2 Kl R 2 ) - ^\2k 1 R 2 )3[{k 11 R 2 ) 
ji(K Il R 2 )fl 1)) '(2K I R 2 ) - ^Ji 1 \2^R 2 ) 3 [{n ll R 2 ) 



«n ?(2), 



(65) 

(66) 
(67) 
(68) 
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With (jHSJ, §23) or (JH7)) . (jHBl) the amplitudes of the 2 linearly independent solutions in region I 
are represented as a function of the energy E. The eigenvalues can now be obtained from the 
ratio of the two continuity conditions (|49|) and (|48|). i.e. from the equation 



dr 



r=Ri 



'd In i? (r) 
dr 



(69) 



This is a transcendental equation for the eigenenergies. It has infinitely many solutions close 
to the energy E = due to the "infinite range" of the Coulomb potential. In spite of the 
complicated nature of the wavefunctions it is not difficult to find the numerical solutions of 
). More explicitly, the logarithmic derivatives in ()69|) look as follows: 



~d\nR (r)~ 






dr 


r= 


Ri 




'dlnRji 


r) 






dr 




r=Hi 




'dlnRji 


r) 






dr 




r=Hi 



^((2Ar)'t/(a ,2Z + 2,2Ar)e- Ar ) 
(2Ar)^(a ,2/ + 2,2Ar)e~ Ar 



'|(^./ i (1) « + S I /f»(r)) 



--Ri 




(70) 



(71) 



=-Ri 



(72) 

Eq. (f?2*|) and (fTTf deal with the case of the exact potential and of the approximate constant 
potential in region I. In calculating the derivative '(r) or /j '(r) one can use properties of 
the different functions involved. 

Having determined the eigenvalues from Eq. (}6T?j) one obtains the amplitude Co from Eq. (|4"8j) 
and thereafter the correct normalization of the entire function must be achieved by numerical 
calculation of the normalization constant (|52|). 



3 Polarization effects in a macroscopic description 

In Sec. 2, only the polarization of the insulator which is produced by a constant positive surface 
charge density Ze$/ {AirR^) on the surface S 2 of the metal was taken into account. However, 
the charge distribution p e (r) of the externally localized electrons produces not only additional 
polarization effects in the insulating surface layer but also modifies the distribution of the 
positive surplus charge on the surface S 2 of the metal. In this section, we present a simple 
macroscopic description of the polarization effects which is based on Maxwell's equations in the 
presence of homogeneous dielectric media. 

Outside of the metallic core, i.e. in the region I (R 2 < r < Ri) of the insulator, and in the 
vacuum (r > Ri) surrounding the insulator the "electric field" E and the "displacement" D 
are related by: 

D(r) = eiE{r) = (1 + Aitxi)E{r) for R 2 < r < R 1 region I , (73) 
D(r) = E(r) for r > Ri region , (74) 

where e\ and \i are the "dielectric constant" and the "electric susceptibility" of the insulator, 
respectively. The fields satisfy Maxwell's equations 

div D = 47rp e (f) , (75) 
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rot E = . (76) 

For the sake of simplicity, we assume that the distribution p e (r) of the externally located 
electrons and, as a consequence, also the positive surface charge distribution ps 2 if) 011 ^2, 
are azimuthally symmetric, i.e. that they depend only on the polar angle •& and the radial 
coordinate r. They can thus be written in the form 

Pe(rO=PeM) = £Mr)*«>(#) > ( 77 ) 

0S2 



1=0 

00 

PS2 (f) = ps 2 (R 2 ,#) = £pf 2 • (7E 

z=o 



Integrating the surface charge density over the surface S 2 must yield the total positive surface 
charge Ze$. This implies the condition 

Ze = J d 3 r 5(r - R 2 ) p s , 2 (R 2 , 0) (79) 

or, with (JZHJ) 



P?=o = PP^ ■ ( 80 ) 



The components pi(r) of the electronic charge distribution are related to the total electronic 
charge distribution by 

pi{r) =2tt J d#sin$p e (r,$)Y l0 ($). (81) 


The Maxwell equations (J7H|) and (j^Sj) are equivalent to the Poisson equation 

A$i(r,0) =-— p e (r,t?) (82) 
£1 

in region I and to 

A$ (r,tf) = -4vrp e (r,$) (83) 

in region 0. Here, the potential 

$(r, t?) = $j(r, t?) for i? 2 < r < R x , (84) 
$(r, 1?) = $ (r, 1?) for r > i?i (85) 

is related to the electric field E by 

E = -V$ . (86) 
The following boundary conditions must be fulfilled at the surfaces S 2 and Si: 

,88) 

£l OR, = OR, m 
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= $o(#i,tf). (90) 
Furthermore, the potential $o must satisfy the asymptotic condition 

$ o ( r>1? )_>0(T) . (91) 



The conditions (|87jl and ()89jl are obtained from the Maxwell equation (|75jl . and the conditions 
(JKKJ) and (|5Dj) result from (|75|) and the requirement that the classical electric field should be 
everywhere finite. 

The solutions $i and $o of the Poisson equations (}82|) and (}83j) in the regions I and 0, 
respectively, can be written in the form: 

$ I (r,i?) = — tf(r,tf) +(pi(r,#) , (92) 
$oM) = tf(r,0) + p o M), (93) 
where \I/(r, $) is the Coulomb potential produced by the electronic charge distribution p e {f ^ 

*W)'-=fMy^' ( 94 ) 

The charge distribution p e (f) itself is to be calculated from the wavefunctions of the bound 
and predominantly externally located (i.e. in region 0, r > Ri) electrons. The "conduction 
electrons" which are located internally (i.e. in region II, r < R?) do not contribute to p e (r) 
because they are already taken into account by the number Z of positive surplus charges. The 
density p e (f) is thus mainly located in region and it is very small in region I, i.e. inside of the 
insulating surface layer, because the eigenenergies of the electrons are far below the positive 
barrier in region I. 

The function ipi(r, $) is the general azimuthally symmetric solution of the Laplace equation 
in region I 

oo / rjl \ 

^iM)=£ C^ + 3TT )M#) (95) 
i=o \ r ) 

and ifo (r, i?) the corresponding most general solution of the Laplace equation in region satis- 
fying in addition the asymptotic condition (|9~Tj): 

OO T-\0 

p»M) = Ei^). (96) 

1=0 ' 

The Coulomb potential \I/(r, i?) can be written as 

00 Ait ( , b,(r)\ 
*(r, *) = E [ a ^y + T^T j > (97) 

where the functions a^(r) and &;(r) are related to the components p/(r) of the electronic charge 
distribution by 



ai (r) := J dr' Pl (r')r n - 1 , (98) 

r 
r 

6,( r ) : = / rfr'pz(r')r' +2 . (99) 



o 
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Substituting the explicit forms (jHSJ) to (fHTj) into the boundary conditions (f8*Tj) . (jEHJ), and (fHUJ) . 

we obtain the relations 



4 * ^(W- 6,(iW(i + i; 



pi+2 



+ —^1^ = 



z=o 



«+2 



YM = o 



Z=0 



(100) 
(101) 
(102) 



Furthermore, using the identity (see for instance Ref. for Z > 1 



sin-$ 



dY, 



10 



1(1 + 1) 

V2TTT 



Yi+i,c 



Yi-i, 



(103) 



together with flH2J), flU and flEJ), the condition 
relation 



V2T+3 V2T^T. 

can be shown to be equivalent to the 



4tt 



(21 + l)e r 



ai(R 2 )R l 2 + 



k(R 2 ) 
Ri +1 



+ C}Rl+ D] 



r>l+l 







(104) 



where / > 1. Eq. (|1U4|) is seen to imply that the potential $i(r, •&) is constant on the surface 
S 2 - For I = 0, we obtain from (jl00|) 



Dl = ~[-b (R 2 ) + Rlp s >] 
£1 



and from (jlOlj) and ()102|) the relations 



• 



Hence, with ()105j) . we find 



and 



C l + — = — 

R\ R\ 



Cl = ^-^[-bo(R2) + Rlp^} 



D° = An[~b (R 2 ) + R 2 2 p s ^} 
For I > I, we obtain from (fffinj) and (jTTHjl : 



a 1 



A 1 



47T 



(2/ + l) £l 

47T 

(2/ + l) £l 



S2 



rt 2 



For the coefficients Df, we find the relations 



(105) 

(106) 
(107) 

(108) 
(109) 

(110) 
(111) 

(112) 
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and 



1 (Z + l) 



(113) 



from the Eqs. (ITU2J) and (fTHTl) . which hold for I > 0. 

The equality of the righthand sides of ()112|) and ()113|) together with (jllOj) and ()111|) yields 
a condition for the surface charge components pf 2 for Z > 1. From this equation one finds the 
explicit form 



pf 2 = ^~ X ~ - ~tt^ ' r '\;;/ r— - (114) 





) + k(R 2 )n: m+1 > ■ < 




1 




+ (*) 


2Z+1 "I 


\ 



for Z > 1. The component p 2 is determined by the conservation of the total positive charge on 
S2 and is given by Eq. (JHU|) . 

For Z > 1, one obtains from (fTTUJ) . JTTTft and (fTT3l 



As we have mentioned already, the charge distribution p e (r, 1?) of the externally located electrons 
is vanishingly small for r < R\. Consequently, we have 

bi(r) ^ for r < R 1 (116) 

and 

a;(r) « a,(i?i) for r < R l . (117) 

The potentials $0(^)5 ^11 — ^1(^2), and the charge distribution p S2 (R 2 ,$) on the 

surface S 2 are thus fully determined for a given azimuthally symmetric distribution p e (r,$) of 
the externally located electrons. 

The potential V(r) acting on an electron in our "mesoscopic atom" is related to the potential 
functions as follows: 

V(r, 0) = -e $ n - V" = -e ^(R 2 ) - V* for r < R 2 , (118) 
V(r, 0) = -e $i(r, 0) + 5 for i? 2 < r < R l , (119) 

V"(r,i?) = -e <&o(?", 0) for r > R 1 (region 0) . (120) 

The potential (— V^) represents the interaction of a conduction electron with the average ionic 
background charge in the metal and has to be added to the constant potential — e $n produced 
by the positive surplus charge on S 2 - The potential barrier B to be added in the insulating 
surface layer is related to the average energy gap between the electronic bands. Since it is of 
quantum mechanical origin, it is not contained in the classically determined term — eo$i(r, 1?). 

The Schrodinger equation with the potential V(r, 0) can no longer be solved analytically. 
One could diagonalize the Hamiltonian 

H = -^A + V(r,$) (121) 

within the set of analytical eigenfunctions which we have determined in Sec. 2. This represents 
a selfconsistency problem, because the eigenfunctions of the Hamiltonian (j!21|) depend on the 
charge distribution p e (r) which is calculated from them. The solution thus requires an iteration 
process. 

Let us consider two special cases for which the results of this section simplify considerably: 



15 



1. If we neglect the polarization effects produced by the external electrons, that is if we put 

p e (r) = 0, (122) 

the quantities ai(r), h(r) vanish for I > 0, and, as a consequence, so do the components 
1 . The positive surplus charge density on S 2 is then spherically symmetric and given 

by" 

P5 2 (r,^)=p 52 Foo = ^||. (123) 

It is easily seen that, in this case, the potential V(r, $) becomes equal to the potential 
V(r) we considered in Sec. 2. 

2. If we assume that the charge density p e {r) of the external electrons is spherically sym- 
metric 

p e (r)=Pe(r), (124) 
only the component pi=o(r) in (fT7|) is unequal 



Pi{r) = S l0 ■ V47rp e (r) (125) 

and, consequently, the quantities a/(r), &z(r), and pf 2 vanish for I > 1. In this case, the 
potentials $n ; $1 and $0 take the following form 

$11 = $i(i? 2 ) (for r < R 2 ) , (126) 

*i(r) = (a (r) + ^ W + + (for R 2 < r < R 1 ) , (127) 

V r J ei Rt sir 

$ (r) = [ a (r) + ^ ] v 7 ^ + — (for r > R{\ . (128) 
\ r J r 

The case of a spherically symmetric distribution of the externally located electrons is realistic 
at finite temperature, because magnetic substates of a given electronic level are expected to be 
occupied with the same probability whenever the mesoscopic atom is in contact with a thermal 
reservoir of finite temperature. 

The functions ao(r) and b (r) depend themselves on the electronic distribution. For r 3> Ri, 
a (r) decreases to zero as one can see from (|UH|) . The quantity \/47r6 (r) represents the total 
charge of the externally located electrons inside the sphere of radius r 

r 

Q e (r) := / dr'r' 2 J dQ' p e (r' \ . (129) 


Using (|ZZ|) in (fT^jl we find 

V^b {r) = Q e {r) . (130) 

For r >> R\, ^/4irbo(r) tends to the total charge Q e (oo) of the externally located electrons. 

If the mesoscopic atom is neutral, i.e. if there are as many externally bound electrons as 
there are positive surplus charges on the surface S2 of the metallic core, the term b (r)\/47i/r 
in Eq. (j!28J) cancels the potential Ze /r for r 3> Ri, as it should. 

Let us study the asymptotic behaviour of the function ao(r) and bo(r) for r — > 00: For 
simplicity, let us assume that only the lowest external electron state is occupied by two electrons, 
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one with intrinsic spin up and the other one with spin down. The external state of lowest energy 
belongs necessarily to orbital angular momentum (I = 0) and to zero number of radial modes. 
We assume that the wavefunctions we determined in Sec. 2 are realistic for r 3> R\. In region 
0, the electronic density distribution is then given by 

Mr) = (-e„) ■ 2RUr; E) i£ = (^g^f (131) 

or, substituting the explicit wavefunction (J2DJ), by 

Pe (r) = t^ClM^ 2, 2Ar) . (132) 

Here, M is the overall normalization constant introduced by Eq. ([52)1 . 

For 2 ^> 1, the function U(ao, bo, Zq) can be replaced by (see Ref. [E], Eq. (13.5.2)) 

TT( h \ ~ 1 f -I ■ ^ (Qo)n(l+a -&o) w / 1 \»\ 

Retaining only the 1 st term in ()133j) . we obtain for the asymptotic density (2Ar <C 1) 

p e (r) » (- e o)(CoAA) 2 (2Ar)do _ 2e _ 2Ar 

27T 

and for the function a (r) 

ao(r) « V^Jdr'r' Pe (r') = ^^f^J ] dz^z^e^ (135) 
The integral in ()135j) is the "incomplete T-function" (see Ref. Eq. (6.5.3), Re(a) > 0) 



DO 



r(a,x) := / dtt^e- 1 (136) 



with the asymptotic expansion (Ref. [H], Eq. (6.5.32) 



!■(«,«) » x-'e- (i + i^l + - ^ - 2 > + . . 1 . (137) 



X X 2 



From (|135|) - ()137j) . we obtain ao(r) in the form 

«o(r) - ^|^ (2Ar)*- 1 e -^ (138) 

for 2Ar > 1. 

Of course, one could raise the accuracy of the result ()138j) by incorporating one or two 
more terms in ()133|) and (|137|) . The important point, however, is that ao(r) tends to zero 
exponentially for large values of r ^> Ry. This justifies the assumption that the electronic charge 
density p e (r) in the asymptotic regime r 3> Ry is correctly described by the wavefunctions of 
Sec. 2, which do not contain the electronic polarization effects. 
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The function bo (r) can be written in the form (see (JSHJl) 



b (r) = J dr'r' 2 p e (r') - J dr'r' 2 p e (r') , (139) 

r 

with 

oo oo . . 

/ dr'r' 2 p (r') = [ dr'r' 2 p e (r') = ■ ( 140 ) 

J •/ V47T 



The 2 nd term one the r.h.s. of ()139j) is easily seen to have the approximate form 

dr'r' 2 Po (r') « ( ~ e ^°^ )2 (2Ar)^ e ~^ (141) 
for 2Ar > 1. 

Thus the function &o( r ) is asymptotically (2Ar > 1) given by 

6 »<^w + ^w (2Ar)d "'^- (142) 



which means that it approaches exponentially the asymptotic value Q e (oo) / '\/47r. 

In the general case of a ^-dependent electronic charge density p e (r, the quantities a^o( r ) 
and bi^o(r) are different from zero and their asymptotic behaviour can be investigated in an 
analogous way. One finds that ai >0 (r) tends exponentially to zero and h >0 (r) exponentially 
to constant values as r goes to infinity. Therefore, at large values of r, the dominant part of 
the potential V(r, i?) is the potential V(r) given by Eq. (jSJ). Thereby, our assumption that the 
wavefunctions determined in Sec. 2 are essentially correct at large values of Ar is seen to be 
justified. 

For non-asymptotic values of r, the electronic charge density and the resulting functions 
ao(r) and bo(r) must be obtained from the eigenfunctions of the Hamiltonian ()120j) which 
contains the potential V(r, i?) defined in ()118|) to ()12()|) . Let us formulate the special case of a 
spherically symmetric charge density and the resulting spherically symmetric potential V(r) in 
somewhat more detail: 

The electronic charge density at a finite temperature T is given by 



/ 



n=i 1=0 [ exp ( + 1 1 



m=— I 



or, using the relation (see for instance Ref. p 7 ], Eq. (3.69)) 



E |^ m (^)| 2 = ^^ (144) 



m=— / 

by 



Pe(r) = -co E E ■ 7 — n ■ E ^ • ( 145 ) 

n=u=o 27r {exp(^^)+l} 

Here, _E n ; is the n-th eigenvalue for given angular momentum I and R\{r\ E n i) is the correspond- 
ing eigenfunction of the radial Schrodinger equation. The parameter p represents the chemical 
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potential which corresponds to the number of externally located electrons. 
The potential V(r, i?) has the form 

V{r) = -eo$i(i? 2 ) - for r 2 < R 2 , (146) 

V-(r) = -eo^+Bo 

= B Q -eA— (a (r) + ^]+C l + ^)Y 00 for R 2 < r < R x , (147) 



Ei \ r 



V(r) = -e $ (r) = -e |4vr^ao(r) + ^M^ +^|f o for r > R 1 , (148) 

where D\, Cq, and are given by the Eq. (jl()5|) . (jl()8j) . and (jl()9|) . resp. The functions ao(r) 
and bo(r) are defined by (fHSJ) and (jHIIJ), where po(r') = \f^ixp e {r'). The wavefunctions Ri(r; E n {) 
are solutions of the radial Schrodinger equation with V(r) given by (jl46|) - (j!48J) 

{"2M (*3 + ;* - V 1 ) + K(r) - E *\ R ^ E ^ ■ < 149 > 

Since V(r) depends itself on the functions Ri(r; E nt ), one has to solve ()149|) by iteration. One 
way to find the solutions Ri(r; E nL ) of (|149|) is to expand the solutions Ri(r; E n i) in terms of the 

o 

radial eigenfunctions Ri(r;Eo) defined in Sec. 2, where the potential V(r) is given by Eq. (|5| 
and related to V(r) by 

V(r) = lim V(r) (150) 

p e (r)-^0 

^(r;^) = Y l CuRi(r;E„) (151) 

V 

o 

.Ev are the eigenvalues obtained in Sec. 2. In Eq. ()151|) we left away the basis functions Ri(r, E) 

o ^^^^ 

ci >rr< ): M ' : n g to continuous eigenvalues E of Eq. Q. Substitution of 1)151)) into ()149)) leads 
to a nonlinear system of equations for the expansion coefficients c u which must be solved by 
iteration. 



4 Results 

We present results on the spectrum of eigenenergies, the eigenfunctions, and on the electro- 
magnetic transition rate between externally and internally localized states for a range of values 
of the parameters R±, R 2 , and £/. So far, we have not taken into account the effects presented 
in Sec. 3, i.e. the polarization of the insulating dielectric material by externally localized elec- 
trons. The additional potential produced by this polarization is expected to modify the wave 
functions in the vicinity of the surface Si, perhaps even rather substantially. At large distances 
(say a couple of Angstrom outside of the surface Si), the additional polarization potential tends 
rapidly to zero. Consequently, the wavefunctions presented in this article are expected to be 
trustworthy at distances of a couple of Angstrom from the surface Si. 

4.1 Eigenfunctions and eigenenergies 

We studied only the case that the Coulomb potential — Ze^/{rei) within the insulator (R 2 < 
t < R\) is small compared to the barrier height B which is related to the gap parameter B$ by 
Eq. (plj. This means that we neglected the Coulomb potential in region I. 
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Table 1: The energies ZeQ/(-R 2 £i), Ze^{e\ — l)/(i?i£i), and the barrier height B as a function 
of the radii Ri and R2 and of the dielectric constant e\. Energies are given in meV and radii 
in A. The parameter B is chosen equal to 2500 meV or 4000 meV. 



R2 (A) 


20 


20 


20 


10 


10 


10 


3.5 


3.5 


JX\ [A) 


98 


98 


98 


1 8 


1 8 


1 8 
10 


1 n 

1U 


1 n 


£\ 


115 

(Ti0 2 ) 


28 
(HgO) 


10 

(AI2O3) 


115 

(Ti0 2 ) 


28 
(HgO) 


10 

(M203) 


115 

(Ti0 2 ) 


10 

(AI2O3) 


Z 


1 


1 


1 


1 


1 


1 


1 


1 


S ( meV ) 


6.26 


25.7 


72 


12.5 


51.4 


144 


35.8 


411 


%e, ( meV ) 


509.8 


495.9 


462.9 


793 


771 


720 


1427 


1296 


5 (meV) for B = 4000 


3490 


3504 


3537 


3207 


3229 


3280 


2573 


2704 


B (meV) for £ = 2500 


1990 


2004 


2037 


1707 


1729 


1780 


1073 


1204 



In Table 1 we present the energies Ze^^R^ei), Ze^{e\ — V)j{R\Ei), and the barrier height B 
for typical values of the parameters i?i i2 and e\. It is seen that neglecting the screened Coulomb 
potential Ze^/ (Ri£i) in region I as compared to B is sometimes a very good approximation, 
but, not in general, especially not for Z > 1 and for radii i?i j2 < 10 A. 

The eigenvalues were determined by a numerical solution of the eigenvalue equation (jnHj) 
using the form of d In Ri(Ri)/dRi as given by the Eq. (J7TJ), and Up to the common 
normalization factor M (see ()52)l ). the radial eigenfunctions R(r;l,E) are then obtained from 
the equations (fTTj) . (flSjl. and (j2D|) in the regions II, I, and 0, resp. 

In order to simplify the comparison of results which correspond to different choices of the 
parameters i?i, 2 , B, Vjj and Z, it is advantageous to use the energy unit Ze^/Ri, i.e. the value 
of the Coulomb potential at the outer surface Si. The resulting dimensionless quantities are 
denoted by 

(152) 



' Ze 2 jRi 

h: =Z^Rr' (153) 

v « ■= z^k • ( 154 ) 

(Vu + Ze /R 2 ) E max _ R x , > 

Ze I /R 1 Zefi/Ri R 2 

Furthermore, it is convenient to introduce the typical length parameters Ro and Roo (= 
0,264589 A): 

fi 2 _ Rqo 
2Me 2 Z = ~Z 

It what follows we give some technical details concerning the determination of the eigenvalues 
and eigenenergies. 

In region I, due to our approximation, the radial eigenfunction R\ is given by Eq. (JTSj) . i.e. 
as a linear superposition of the modified spherical Bessel functions f\ (x) and fi 2 \x) defined 



^0 := ^7377 = ■ (156) 
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in (|5H|) and (|54j) with x = Kir. We calculate these functions by representing them in the form 
(see Ref. 0, Chapt. 10) 



n V ' 2x 



R[l+~, -x) e x + {-I) 1 RII + -, x) e~ x 



2x' 



R ( I + —, x ) e 



(157) 
(158) 



as a linear superposition of exponentially rising and decreasing terms. Here, the function 
R(l + §,x) is defined (see Ref. 0, Eq. (10.2.11)) as 



(159) 



Z + -, A; : 



fc=0 

(Z + fc)! 
fc!(Z-Jfc)! 



It turned out to be convenient to rewrite R\(r) in the form 

1 



RJr) = A' l gl 1 \x)e x ~^ R2 + -\B , 1 + {-l) l A > l e- 2 ^ R2 

7T L 

■ g {2 \x)e- { - x - RlR2) , 

where the functions gf~ ,2 \x) are related to fi\x) by 

^ {1) (x) :=/«(x)e-% 
^ (2) (x) :=// 2) (x)e^. 

The amplitudes A[, B{ are connected with the amplitudes A\ and B\ of Eq. (JTSJ) by 

A[ = A ie KlR2 , 

B[ = [(-l) l+1 A I + irBj}e- KlR2 . 



(160) 



(161) 
(162) 



(163) 
(164) 



______ 

The advantage of using ()160|) instead of (f~) is that the functions g t ' (y) can be numerically de- 

(1 2) 

termined more comfortably than the functions /y ' (y) and that they have a simpler asymptotic 
behaviour 



2.T 



7T 



2a; 



From the continuity conditions ([46)1. (}4~l one obtains and -Bj in the form 



4 
3 



7T 



(ACli2 2 



dgi 2) (KiR 2 ) -ji{KiiR 2 ) - —gl 2) {KiR 2 )j'i{KiiR 2 



Ki 



^ g m {KlR2 )f {KllR2 ) 



(165) 
(166) 

(167) 
(168) 
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where the function dg^ 1,2 \x) are defined by 

dg^ix) = e~* ff Y (x) , (169) 
dg®(x) = e**fF>\x) . (170) 

Then, the continuity condition (Jfinj) assumes the explicit form 

F(e) : = U \l + 1 - y , 21 + 2, 2XR^j [IR 1 + k^Rj] 

+ 2Afli-(Z + l-^)E/^ + 2-^,2Z + 3,2Ai^ ■ Ri (171) 
+ (Ai2i-0c(i + l-y,2Z + 2,2Ai2 1 J -i2i = 



with 



Ri:=Ri(R!), (172) 
- l ' A' I gj%{K 1 R 1 )e- K ^ Rl - R ^e- 2K ' R2 . (173) 





In terms of the dimensionless parameters (j!52j) to ()155j) the arguments which appear in the 
transcendental equation (j!71j) take the form 

(174) 
(175) 
(176) 
(177) 

We determined the zeros of the function F(e) by calculating numerically F(e) in small steps 
of e between limits e m ; n and e max . Whenever F(e) changed sign from a given step to the next 
one, we refined the mesh by usually a factor of 100 so as to determine the zero with the desired 
accuracy. We note that, for each given value of the angular momentum I, there is an infinity 
of eigenvalues as e approaches e = from below. This is so as a result of the infinite range 
of the Coulomb potential in region 0. In the domain of small |e|, i.e. of very weakly bound 
electrons, and also for large well depth vn, the function F(e) oscillates between large positive 
and negative values. Nevertheless, for a sufficiently narrow system of mesh points, the function 
F(e) is of reasonable size for the meshpoint closest to a zero. Therefore, we encountered no 
problems of numerical accuracy in determining the eigenvalues. 

Once an eigenvalue is determined, the amplitudes A[, B[ are obtained from ()167|) . ()168|) 
then Ai, Bj from (| 163)1 . ()164j) . The amplitude Cq of the radial wavefunction Ro(r) in region 
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(see Eq. ()34j) is found from the continuity condition ()48)1 and, finally, the overall normalisation 
constant M from ()52|). The integral in Eq. (p)2~|) is done numerically. 

We note that the numerical determination of the eigenvalues and eigenfunctions turned 
out to be extremely rapid. Thus it is conceivable that one uses these eigenfunctions (in prin- 
ciple completed by the scattering states) as basis functions for diagonalizing the additional 
polarization potentials studied in Sec. 3. 

We note in passing that A[ has in general a very small value which is difficult to obtain with 
sufficient accuracy from Eq. ()167|) . It turned out that an accurate value of A[ can be obtained 
by solving the equation (jl71|) for A[ at the value of the eigenvalue which had been determined 
before. 

It is found that there are 3 classes of eigenfunctions: "external" or "class l"-states which 
are predominantly localized in the vacuum outside of the insulator (r > Ri), "internal" or 
"class 2" -states which are mainly localized within the metallic core (r < R 2 ), and "class 3"- 
states which exhibit a non-negligible probability amplitude both in region and region II. For 
a sufficiently large value of B, the probability for finding the electron in region I, i.e. within 
the insulator, turns out to be negligibly small. The spectrum of class 1-states resembles the 
Coulomb spectrum. On the other hand, the eigenvalues of class 2-states, in lowest approxima- 
tion, are given by the eigenvalues of an infinite spherical square well of radius R 2 i.e. by the 
boundary condition ii(KuR 2 ) = 0. 

If K11-R2 3> 1, the spherical Bessel function ji(nuR 2 ) can be approximated by 



3i(khR 2 ) » — 5- cos 



Kult'2 

2 4 



178) 



The zeros of this function together with (jl2j) lead to the energy formula 



„ Ze\ ftV (3 A 2 „„„. 

E '>'-- V »-l£ + 2MRl{i + n+ 2) ' <179) 

where n — 0, 1, ... as well as 

The spacing between neighbouring levels of given I or given n is 

^-^ = 5sJ (§ + *• + a) • W 

E ^- E -' = 2KmA l+n+ V ■ (181) 

One should notice the completely different spacing of the eigenvalues of internally located 
electrons from the Coulomb dominated spacing of externally located electrons. Of course, if 
KuR 2 is not large compared to 1, the simple results (j!78|) to (jl81|) do not hold. 

If it happens that, for given I, the energy of a predominantly externally and of a predomi- 
nantly internally localized state come very close to each other, the corresponding eigenfunctions 
are not negligible in both the regions II and 0. We call eigenstates of this nature "class 3"- 
states. These states play an important role in the lifetime of the system and we shall come 
back to them in Subsec. 4.2. 

In Fig. 3, we show the spectrum of eigenvalues e for I = 0, I = 1, and I = 10 for a given 
set of parameters. We note that, apart from the dimensionless parameters defined in the text, 
lengths are given in Angstrom (10 -10 m) and energies in meV (10~ 3 eV). 
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Figure 3: Comparison of the spectra of I = states (Fig. 3a) and I = 1 states (Fig. 3b) and 
I = 10 states (Fig. 3a) for the following set of parameters: Z — 1, Ri — 28 A; Ri — 20 A; 
Vu = 15950 meV; b = 5.83; v n = 31.031, R = 0.271 A; e min = 0.1; e max = 32.41; energy mesh 
Ae = 0.0001. Energy unit: Zel/R 1 = 514 meV. 



Given the energy unit Ze^/Ri = 514 meV in Fig. 3, the |e|-values belonging to externally 
located states must be smaller than 514 meV. The lowest-lying external state corresponds to 
the lowest line of the bunch of narrowly spaced lines at the top of the spectrum at e — 0, 68 
or E = —350 meV. Note that the energy difference between the first excited external s-state 
and the lowest external s-state is about 69,5 meV. This is also the order of magnitude of the 
average distance between low- lying external levels. Since the "room temperature" of 300° K 
corresponds to 26 meV, one may conclude that, at room temperature, for the parameter set 
chosen in Fig. 3, the externally bound electrons are predominantly in their lowest available 
states. 

The spectrum of eigenvalues looks similarly for different fixed values of the angular momen- 
tum I. The fact that we show only the subsets of eigenvalues for given /-values in Fig. 3 should 
not lead to the erroneous impression that the levels of internally located electrons are so widely 
spaced that in general, no level of internal electrons is to be expected in the range of energies 
of external electrons. This is not so due to the fact that there are also the levels for I ^ 0. 
Consequently, there are generally " intruders" of internal levels in the spectrum of external ones 
and, consequently, there is in general also a limited number of "class 3-states" which, as we 
shall show, lead to rapid decay. 

We should also notice in Fig. 3 that the distance between successive eigenvalues for given 
/ increases with increasing eigenvalue, i.e. as \e\ — > 0. The reason for this trend is qualita- 
tively seen from the formulae (jl80|) and (|181|) . For the parameter set chosen in Fig. 3, these 
approximate formulae are only valid for the high-lying levels. 

In Fig. 4, we show the wavefunctions corresponding to internal s-states with 0, 1, and 8 
radial nodes. It should be noted that even for the case of the highest one with 8 radial nodes 
at |e| = 17, 96, the localization within region II is perfect and the tails in region I are very tiny. 

In Fig. 5, we present a comparison of external s- wavefunctions for radial node numbers 
n — 0,1, and 5. Again one notices that the wavefunctions dive only weakly into region I and 
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Figure 4: Internally located wavefunction (class 2-states) with radial node number n = 
(Fig. 4a), n = 1 (Fig. 4b), and n = 8 (Fig. 4c). Orbital angular momentum I = 0. Energies 
(in units of Zeg/i?i = 514.3 meV): |e(n = 0;/ = 0)| = 32,235; |e(n = 1;/ = 0)| = 31,698; 
\e{n = 8; I = 0)| = 17.963. Choice of parameters as in Fig. 3. 




Figure 5: Externally located wavefunctions (class 1-states) with radial node number n = 
(Fig. 5a), n — 1 (Fig. 5b) and n = 5 (Fig. 5c). Orbital angular momentum / = energies 
(in units of Ze\jR x = 514.3 meV): \e(n = 0;l = 0)| = 0.63325; \e(n = l;l = 0)| = 0.44533; 
|e(n = 5; / = 0)| = 0.18213. Choice of parameters as in Fig. 3. 
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are completely negligible in region II. The energy difference (E n= i t i =Q — E n=0t i =0 ) between the 
two lowest-lying external levels turns out to be 96 meV. The parameter set was chosen to be 
the same for the Figs. 3-5. 

For the high-lying eigenenergy E n= $ ^ = o ~ —93, 5 meV the most external maximum of the 
wave function is seen to occur at r ~ 130 A whereas the wavefunction of the lowest external 
s-state with n = becomes vanishingly small already at r ~ 60A. This strong dependence of 
the radial extension on the energy of the state is typical for a Coulomb-spectrum. 

The class 1- and class 2-states are thus of very different nature. The eigenvalues of class 
1-states must be above the potential energy (— Ze^/Ri) on the surface Si of the insulator. 
Their spacing is much smaller than the one of class 2-states and tends to zero as the binding 
energy approaches zero. The lowering of the well-depth —V^ by the potential produced by 
the positive surface charge on the metal surface S 2 and the screening effects in the insulator 
is given by — \Ze\j (eiR 2 ) + Ze^(ei — l)/(eii? 1 )]. This lowering of the Fermi energy E F in the 
metal is larger than the minimal binding energy — [Ze^/ Ri] of the lowest external electron by 
the amount 
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For the choice R\ = 28 A and R 2 = 20 A and the charge number Z = 1, we have 
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i.e. a number between -20 and -200 meV depending on the dielectric constant e\ of the insulator. 

If there are unoccupied levels between the Fermi energy Ep and (Ep + |Ai?|), an electron 
in the lowest external state can decay into a lower-lying internal state. 



4.2 Electromagnetic decay of external states 

For the sake of simplicity, let us consider a system with one external electron and Z positive 
surplus charges. If the external electron is in some arbitrary externally located state it will first 
decay to the lowest-lying external state by a cascade of unhindered transitions. This lowest- 
lying external state is characterized by the angular momentum I = and has no radial nodes 
(n = 0). 

Subsequently, this lowest class 1-state will decay into lower-lying unoccupied internally 
located states (class 2-states). The strongest transitions are of multipolarity El and we only 
consider such dipole transitions. Even these most rapid transitions between a class 1 and a 
class 2 state are reduced in comparison with transitions between states of the same class due 
to the small overlap between externally and internally located wavef unctions. 

The probability per time-unit for a dipole transition from the initial state tpi to the final 
state ijjf is given by |Hj 

^^tlfcT 1 ^ 1 ^ 12 ' (183) 

where E2 and E\ denote the final and initial energy, resp., and a = e\j (he) is the fine structure 
constant. 

The initial and final wavefunctions are 

= R h (r,E 1 )Y hmi (^<p) , (184) 
Mr,#<p) = R l2 {r,E 2 )Y hm2 ('d, ( p) , (185) 
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Table 2: The quantities W z and W± entering the dipole transition probability (see Eqs. (J188j) . 



h-h 



h-h 



m 2 = ni\ + 1 



W. 



(?i-mi-l)(?i-mi) 
(2h+l)(2h-l)-2 



(h+mi+l)Gi+mi+2) l 2 
(2h+l)(2fi+3)-2 



7712 = mi 



W z 



{h-mi)(h+mi) 
(2fi+l)(2Zi-l) 



w 2 



(fi-mi+l)(fi-mi+2) l 2 
(2fi+l)(2h+3) ~ 



m 2 = mi 



(;i+mi)(^i+mi+l) 
(2fi+l)(2h-l)-2 



{h-m 1 +l)(h-m l +2) 
(2fi+l)(2h+3)-2 



where we leave away the spin-functions, and where , i?; 2 are the properly normalized radial 
functions 1)51)1 . The matrix-element factorizes 



(^|fIVi>= drR l2 (r)r 3 R h (r)- dQY*e r Y hmi 



(186) 



It is easily seen that the following relation holds: 



y rff2Y" iimi e r Y/ 2J)12 



\w + \ 2 + \w_\ 2 + \w z \ 2 , 
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where the quantities W± and VF Z are defined by 



W ± := T^f J dQY* mi Y h±1 Y hm2 
W K :=J^ [ dnY* mi Y 10 Y hm2 . 



(188) 
(189) 

They are evaluated in terms of Clebsch-Gordan coefficients using the Wigner-Eckart theorem 



W ± := Tl 



I2l 2 + 1 
2Zi + l 



C(/ 2 l/i, m 2 , ±1, mi)C(Z 2 lii000) 



W z : = y|^C(/ 2 l/ 1 m 2 0m 1 )C(/ 2 l/ 1 000) . 



(190) 
(191) 



The non-vanishing values for W± and W z are presented in Table 2. The radial integral, i.e. the 
1 st factor in Eq. (J186j) is evaluated numerically. 

We consider the case that the initial state is the lowest external state, i.e. a nodeless 
wavefunction with l\ = 0. Thus, the final angular momentum / 2 is equal to 1. 

The final magnetic moments m 2 may have the value m 2 = —1, 0, +1. 

We are mainly interested in the lifetime r of a "mesoscopic atom" which is the inverse of 
the total decay rate P: 

(192) 



t:= p- 



The total decay rate P is defined by summing the transition probabilities Wu over all the 
unoccupied final states which are connected with the initial state by a finite matrix 
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Table 3: Separation energy Wq of an electron, Fermi-energy E F and well depth of uncharged 
metals (in meV). 
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Table 4: Lifetime r as a function of the outer radius Ri (A) for: Z = 1; R 2 = 20 A; B = 
3000 meV; V$ = 15950 meV (Al). Lifetimes in units: sec, minute (min), day (d), year (a). 
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element (^1 1*^1 ^2) - 

p = Y^'W 12 . (193) 
2 

Considering only dipole transitions is, of course, an approximation which yields only a lower 
limit for the total transition rate and an upper limit for the lifetime. 

The restriction of the sum in ()193|) to unoccupied internal electronic states is indicated by 
the prime '. 

At temperature 0, a dipole transition starting from the initial external s-state can go into 
all the l 2 = 1 states the energy of which is between the Fermi energy E% and the energy E\ of 
the initial state 

E F < E 2 < Ei . 

At finite temperature, one has to weigh all the final states with the factor 



exp(^p)+l_ 

So far, we calculated the lifetime for the case T = only. 

In Table 3, we present the well depth V^, the separation energy W , and the Fermi energy 
E F = Vjj — Wq for a couple of metals. We note that these empirical parameters have been 
determined for materials with macroscopic size, not for mesoscopic clusters of such materials. 
Nevertheless, the parameters can be expected to be essentially correct also for mesoscopic 
systems. 

Let us now come to the results we obtained for the lifetime as a function of the different 
parameters of the system. We studied separately the dependence of the lifetime r on the 
parameters R\, B, Z, and R 2 , keeping always all but one parameter fixed. 

Dependence of r on Ri 

Augmenting R\ for fixed value of R 2 enlarges the width of the barrier and, consequently, 
decreases the amplitude which an externally located electron has in the inner region II. As 
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Figure 6: Lifetime in sec and logarithmic scale as a function of radius R\. Apart from R\ , all 
the other parameters are chosen as in Fig. 3. 
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Figure 7: Lifetime r in sec and logarithmic scale as a function of the barrier height B. Apart 
from B , all the other parameters are chosen as in Fig. 3. 



the decay of a class 1-state into a class-2 state depends strongly on this tail, the transition 
probability must decrease and the lifetime must increase as a function of rising R\. The 
results presented in Table 4 support this expectation. The lifetime r is found to rise roughly 
exponentially with (Ri — R 2 ). This is so because the amplitude of the tail in region II depends 
exponentially on — R 2 ) as one can see from WKB. The dependence of r on Ri is also shown 
in Fig. 6 in logarithmic scale. 

Dependence of r on the barrier height B 

In Fig. 7, we show the dependence of the lifetime r on the barrier height B in logarithmic 
scale. The lifetime increases as a function of rising B. The increase is somewhat smaller than 
exponential as can be expected from WKB. 

Dependence of r on R 2 and on Z 

If one changes R 2 with all the other parameters fixed, one modifies predominantly the 
internally located class 2 states and the corresponding energies. As one can see from Fig. 8, 
the lifetime r as a function of rising R 2 decreases on the average as one would expect from the 
fact that the thickness (R\ — R 2 ) of the surface layer decreases. However, beside the average 
decrease of the lifetime, one observes a kind of resonance structure: Within a change of R 2 by 
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Figure 8: Lifetime r in sec and logarithmic scale as a function of R 2 . Apart from R2 , all the 
other parameters are chosen as in Fig. 3. 
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Figure 9: Probability (W) of lowest external electron to be found in the inner region II (R < R2) 
as a function of R2. Apart from R 2 , all the other parameters are chosen as in Fig. 3. 



1 A, the lifetime r can change by 10 orders of magnitude. 
How can we understand these huge fluctuations of the lifetime? 

In Fig. 9, we show the probability of an externally located electron to be found in the 
inner region II. This localization probability in region II, which is of crucial importance for 
the transition rate, is seen to show pronounced peaks at precisely those values of R2, where 
the lifetime has very small values as one can see by comparing Fig. 8 and Fig. 9. The large 
variations of the transition rate as a function of i? 2 ( an d of Z) can be simply explained as 
follows: 

As one changes R2 (or Z), the eigenvalues of class 2-states move whereas the ones of class 
1-states remain essentially the same. More specifically, the eigenvalues E n i = —\E n i\ of class-2 
states are lowered, if the inner radius i? 2 or the positive surplus charge Z is increased. Thereby, 
it happens regularly that a class 2-state which has the same angular momentum I as the 
decaying class 1-state, crosses the energy of this class 1-state. If one augments R 2 or Z on a 
sufficiently long scale, such crossings occur several times. Close to the crossing, the external 
class 1-state develops a non-negligible amplitude within region II and the class 2-state develops 
an appreciable amplitude in the external region 0. We call such mixed states "class 3-states". 
If the decaying class 1-state happens to have the nature of a class 3-state, its decay probability 
to lower-lying class 2-states is considerably enhanced as compared to the ordinary situation 
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Figure 10: Lifetime r (sec, log.scale) as a function of radius R\ and R2 (in A) in a 3- dimensional 
plot. Except for Ri and R2, all the parameters are chosen as in Fig. 3. 



where no class 2-state of the same angular momentum is available. We note in passing that 
this enhancement of the transition probability to lower-lying class 2-states does not only occur 
for dipole transitions but also for transitions of higher multipolarity. So far we have not yet 
calculated the transition rates of higher multipolarity, because they would only matter if the 
decaying external state is a class 3-state. 

The phenomenon can be considered as an example of the mechanism which Landau and 
Zener had studied in famous papers |10j . This is demonstrated in Appendix. 

In Fig. 10, the dependence of the lifetime r on the radii R\ and R2 is shown in a 3- 
dimensional plot. Again one notices that the dependence of r on R\ is smooth whereas the one 
on R 2 contains "resonances". 

In Fig. 11, the energy difference between the lowest-lying class 1 state and the closest lying 
class 2-state of the same angular momentum is shown. The plot supports the explanation of 
the resonance phenomenon in terms of a coupling between close-lying states. 

In the numerical results presented in this section, the values B of the barrier in region I 
and the well depth Vu of the average potential in region II and the dielectric constant e\ were 
chosen arbitrarily, and not determined for given choices of the material. 

5 Summary and Discussion 

In this paper, we determined the quantum states and eigenenergies of an electron in the average 
potential produced by a positively charged metal cluster which is covered by an insulating 
surface layer. Because of the similarity of the externally localized electronic states with the 
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Figure 11: Energy difference between lowest-lying class 1-state and the closest lying class 2- 
state of the same angular momentum as a function of R 2 (A). Energies are given in units of 
Ze^/Ri = 514.3 meV. Except for R 2 , all the other parameters are chosen as in Fig. 3. 



bound states of an electron in the Coulomb potential of an atom, we occasionally referred to 
the system as a "mesoscopic atom" or, in case of a non-zero total charge, as a "mesoscopic 
ion". 

In particular, we investigated the electromagnetic decay of an externally located electron 
("class 1-state") to a lower-lying unoccupied electronic state located within the metal ("class 
2-state"). It turned out that the life-time depends sensitively on the radius R 2 of the metal 
core and on the number Z of positive surplus charges of the metal. In the model we studied, 
we found life-times varying between 10~ 10 s and years depending on the value of R2 or of Z. 

The value B of the barrier height was 3 eV in our calculations. The lifetime depends strongly 
on the value of this parameter and decreases for decreasing B. 

The origin of the large fluctuations was found to be the appearance of states with non- 
negligible amplitudes both outside and inside of the insulating surface layer ("class 3-states"). 
The mechanism leading to these "class 3-states" is described in Appendix. 

The wavefunctions and eigenenergies were calculated neglecting the polarization of the insu- 
lating surface layer by externally located electrons. The potential produced by these electronic 
polarization effects was determined in Sec. 3. It was shown that this additional potential mat- 
ters mainly in the vicinity of the outer surface Si of the insulator and decreases exponentially 
outside of this surface. The polarization potential could be taken into account by diagonalizing 
the total Hamiltonian including the polarization potential in the basis of the wavefunctions 
determined in this paper. 

Apart from this physical simplification of the problem we neglected the Coulomb potential 
produced by the positive surplus charge within the insulator. This approximation is justified 
whenever the change of the Coulomb potential within the insulating surface layer is much 
smaller than the barrier height B . This condition was fulfilled in all the cases we studied 
numerically. 

Let us now turn to the discussion of some open questions: 
Is the model sufficiently realistic? 
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Certainly, our model should only be regarded as a first step which has to be improved in various 
respects. 

The 1 st improvement we already mentioned is to incorporate the additional potential be- 
tween the external electron and the polarization cloud it produces in the insulator. The form 
of this electronic polarization potential was derived in Sec. 3. 

A 2 nd improvement would be to replace the constant potentials — Vq in the metal and +Bq 
in the insulator by periodic potentials defined by the crystal structure of these materials. Of 
course, this is much more complicated and, to our knowledge, it has rarely been done so far for 
mesoscopic systems. Furthermore, the crystal structure in mesoscopic aggregates is in general 
poorly known. 

As a result, the electronic eigenvalues of the total system are expected to exhibit a band 
structure similarly to a macroscopic solid, the insulating surface layer being characterized by a 
completely filled valence band and an empty conduction band. The empty conduction band in 
the insulator can be at negative or at positive energy if the energy of the neutral system plus 
one electron at infinity is put equal to 0. Describing the average potential in the insulator by a 
positive barrier Bq > is only meaningful if the conduction band of the insulator is at positive 
energy. 

As a further restrictive condition for our model to be a reasonable 1 st approximation there 
should be no localized traps for an additional electron in the uncharged insulator. This implies 
that the individual atoms or molecules of the insulator should not have bound states with an 
additional electron as a monomer [H]. In general, bound states of an electron with a neutral 
atom or molecule do not exist for atoms or molecules with completely filled shells. Therefore, 
one should use insulators which exhibit this property. Presumably, Al 2 03 fulfills this condition. 
Finally, we neglected effects of the electron-phonon coupling altogether. We presume that they 
tend to reduce the lifetime, but perhaps not drastically. 

What were the physical motivations of this investigation 

As long as one can only produce a limited number of "mesoscopic atoms", one could just 
investigate experimentally the effects which we studied in the present paper, especially the 
amusing strong dependence of the lifetime on the radius and the charge of the metallic core. 
Although we would be highly pleased by an experimental scrutiny of our predictions, still more 
exciting open questions could be asked, if it were possible to produce longlived systems of the 
kind we presented in this paper, and in large quantities. In this case, one would have access to 
the study of molecules formed from mesoscopic atoms and to condensed matter consisting of 
mesoscopic atoms as building blocks. 

The presence of externally localized electrons would represent a fundamental difference from 
the case of a condensed system of uncharged mesoscopic clusters. On the one hand, one may 
expect that, in a condensed matter phase of mesoscopic atoms, the weakly bound external 
electrons would partly dissociate from their mother atoms and move freely in the average field 
produced by the charged insulated cores. We would thus surmise that such a system would be 
a (good) conductor. On the other hand, the fact that the building blocks of such a condensed 
matter carry charge, will have a profound influence on the equilibrium state of the material. 
Probably, its density is lower than for uncharged constituents, where the cohesive forces are of 
van der Waals type and, consequently, of much shorter range. 

How far are such dreams from reality? 
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There are evidently two fundamental obstacles to overcome, a theoretical one and an experi- 
mental one. 

The theoretical problem is the stability of the charged and insulated single system, and the 
experimental one is the production of a large number (10 23 ) of such systems. 

As we learnt from this paper, a mesoscopic system of the kind we investigated can only be 
stable in a strict sense, if there is no empty internally located electronic state below the energy 
of the lowest externally located state which is at an energy slightly above —Ze^/R\. 

The positive surplus charge Ze on the surface S 2 of the metal leads to a decrease of the 
well-depth — > Vn as we have shown in Sec. 3: 



If, for the uncharged metal, the Fermi level is at the energy E F — (with V^, E% > 0), it lies 
at the modified energy 



in the charged metal. 

So the condition for strict stability is that there is no level of an internally located electron 
between E F — Vn < and —Ze\jR\. As the separation energy Wq > of an electron from the 
uncharged metal and the Fermi energy E F are related to V^ by 



we can state that strict stability is achieved whenever there is no level between the energies 
— {W — Zel/ei ■ [I/R2 + (£1 — l)/Ri]} and —Ze^/Ri. This means that, metal with small W 
is favourable for stability. 

In an ordinary metal, the separation energy Wo is of the order of a couple of eV (examples: 
W (A1) = 4, 28 eV; W (Ag) = 4, 26 eV). In the case of a radius R 2 = 20 A of the metal core, the 
distance between neighbouring electronic levels is of the order of 100 meV. This means that in 
the case we chose for our numerical investigations, there are many unoccupied levels below the 
energy of the lowest external electron. However, if we consider a system with an inner radius 
i? 2 ~ 4 A, the spacing of internally located electron levels of given / approaches already 1 eV. It 
is therefore conceivable that one might succeed to model a strictly stable ,, mesoscopic atom". 

A further trick to widen the spectrum of internally localized electrons would be to replace 
the metallic sphere with radius R 2 by a metallic bubble of outer radius R 2 and a thickness of 
1 A to 2 A. In this case the levels differing by the number of radial nodes would have a much 
larger spacing. 

Can such a system be made experimentally? 

This is perhaps possible, if one replaces the charged metal layer by a (charged) fullerene 
6*60, which has a radius of about 3,5 A. The separation energy Wo of an electron from graphite 
is known to be W — 3, 8 eV. We surmise that the separation energy from the fullerene is about 
the same. Of course, it would no longer be justified to describe free electrons in the fullerene 
as moving in a simple square well perpendicular to the sphere S 2 . But the spacing of electronic 
levels is expected to be in the order of eV. 

Finally, one is not limited to the spherical geometry. We may also envisage to consider 
charged nanotubes consisting of a cylindrical surface crystal of C-atoms covered by an insulating 




(194) 




(195) 



V° = E° F + Wo, 



(196) 
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surface layer. The charge could in this case be generated by applying an external voltage to the 
ends of the nanotubes. As there is a lot of activity spent on nanotubes presently, we hope that 
one could succeed to produce small macroscopic amounts of insulated charged nanotubes in a 
not too far future. Of course, the problem of the stability of externally localized electrons with 
respect to electromagnetic decay would be similar for nanotubes as for the case of the system 
we investigated. So we hope that our investigation stimulates the interest of the experts in the 
field of mesoscopic physics. 
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Appendix 



Landau-Zener Model for explaining the appearance of class 3-states 

Let us consider the simplified potential 

Ze 2 

V(r) = -Vu6 (R 2 - r) + B6 (R 1 - r)6 (r - R 2 ) °-6 (r - R 1 ) . (197) 

r 

We decompose this potential as follows 

V{r) = V 1 {r) + V 2 {r) , (198) 
ZeT 



-B - 



o 



(199) 



r 

V 2 (r) = -V u e (R 2 ~r) + B6 (r - R 2 ) (200) 
and define the two Hamiltonians 

H m :=f + V m (r) , (201) 
where T is the kinetic energy operator of radial motion 

? = _1L(^ + 1±J1±A\ . (202 ) 
2M \dr 2 r dr r 2 J 

The eigenstates ip\p, of the Hamiltonians (j201|) 

%^\t)=E^^ , (203) 

H 2 <pW(r) = El\W (204) 

form a non-orthogonal set of basis functions. 

We consider an eigenstate ip of the Hamiltonian H with eigenenergy E 

H:=f + V(r) , (205) 

Hi) = Ei/; (206) 

for the special case that the angular momentum I is the same for all the states and that two 
specific eigenenergies E^> and E^ are very close to the energy E in 



Then, we may expect that the expansion of ip in terms of the basis states (p^ and cpy may 
be limited to the two states with the energies E^\ E^f 1 very close to E v \ 

= AM 1] ) + B„\^) . (207) 

In principle, we would have to sum over all eigenstates cpfi and (pffi on the righthand side of 
(|207|) . but of all the expansion coefficients A u , B u only the ones corresponding to the energies 
E)P , E^f 1 closest to E are expected to be large and all the others to be much smaller. 
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Henceforth, we leave away the subscript v for simplifying the notation. It is convenient to 
introduce the dual basis states 0^'', defined by the orthogonality conditions 

{<p®\<p®)=6v , (208) 

where i and j can be 1 or 2. <pV-> and (p 1 ^ can be obtained by diagonalizing the inverse of the 
overlap matrix 

M {ij) = {<p^\ip^) . (209) 

It is easily seen that the Schrodinger Eq. (|2U6|) together with the assumption ()207|) leads to 
the two coupled linear equations 

^ ( S (2) -E + V™ )\B ) ~ U ' (2iUj 

where the matrix elements are defined as follows: 

vty.= (<p®\v 1>2 \<p®) • (2ii) 

Introducing the abbreviations 

m .= E^ + Vf 2 (213) 
the solutions of the eigenvalue equation 

(£« _ £) (^(2) _ £) _ ^12 . y 2 2i = o (214) 

are given by 

E m i ftp) i ^ 

E± = - \ ] ± y(EW - EW)* + AV^V? . (215) 
Writing the two states \ip) in the form 

|V»±>=A t -{b (1) >±A ± | ¥ »«>} , (216) 

where 

A± := ^ , (217) 
one obtains A from the Eq. (|210|) and the amplitudes A± from the normalisation of \if>±). 

(1 2) (1 2) 

The physically interesting quantities are the probabilities W\ ' and Wl ' to find the 
system described by \ip±) in the substate ip^ or ip( 2 \ These probabilities are given by the 
mean-values of the projection operators 

P (li2):= |^))(^)|. (218) 
It is easily seen that the probabilities are given by the following equations 



(i) 



1 + X+M 



12 



*l p <*> = iTCTTaI • < 219 > 



12 ^ 



<^I P <*> = 1 + 2A> + A| ' < 22 °) 
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R 2 
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Figure 12: Schematic plot of the Landau-Zener crossing of an external and an internal energy 
£(1.2) yielding the class 3-state energies E±. 



w! 2) = ^-|p ( 2)i^-; 



1 - A_M 



12 



2A_M 12 + A 2 . ' 
A_(A_-M 12 ) 



1 - 2A_M 12 + \ 2 _ ' 
where M 12 is the overlap of the two basis states 

M 12 = (v^V 2 )) 

and the quantities A± by 

1 



A 



± 



2V 12 



where 

It is seen that for 



_ A £(1,2) ± ^(A£(1.2))2 + W^V 2 

AE := - E^ . 



21 

2 r ' 



(221) 
(222) 

(223) 

(224) 

(225) 
(226) 



|AE| 2 >4|\/ 1 12 -\/ 2 21 | , 

i.e. far away from the crossing of the energies E^ and E^ 2 \ one of the ratios A± is much larger 
than the other. 

For AE > one finds 

s 



A-j- ~ 



A_ 



AE V 2 21 
V^ + AE 



and for AE < the reversed result 



If 1 

A_i_ ~ -— — + 



v; 12 \AE\ ' 



A_ 



T/21 

K 2 

\AE\ 



(227) 
(228) 

(229) 
(230) 
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This shows that a predominantly externally localized state turns into a predominantly internally 
localized one as the energies E^ and E^ cross as a function of R2 or Z. The coupled true 
states ip± are thus seen to change their nature as a result of the crossing of the energies 
E^ of a class 1 and a class 2 state. Close to the crossing point the states ip± are seen to be 
" delocalized" , containing comparable contributions of the two basis states (p^' 2 \ 

This is in essence a Landau-Zener crossing. In Fig. 12 we show a schematic picture of the 
energies E^ 2 \ E + , and E_ as a function of the the parameter R 2 or Z. 
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